Many-body reduced fidelity susceptibility in Lipkin-Meshkov-Glick model 
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We study the reduced fidelity susceptibility \r for an M-body subsystem of an iV-body Lipkin- 
Meshkov-Glick model with r = M/N fixed. The reduced fidelity susceptibility can be viewed as 
the response of subsystem to a certain parameter. In noncritical region, the inner correlation of the 
system is weak, and Xr behaves similar with the global fidelity susceptibility \ g , the ratio r\ — Xr/Xs 
depends on r but not N. However, at the critical point, the inner correlation tends to be divergent, 
then we find Xr approaches Xg with the increasing the N, and r\ = 1 in the thermodynamic limit. 
The analytical predictions are perfect agreement with the numerical results. 
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I. INTRODUCTION 

Quantum phase transition (QPT) [l[ occurs at ab- 
solutely zero temperature is driven purely by quantum 
fluctuations. It was studied conventionally by Landau 
paradigm with order parameter in the frame of statistics 
and condensed matter physics. Recently, two quantum- 
information [pi concepts, entan glem ent El, 4 , H , M, 3, HI , 
1 GS El El El and ndelityjil El, EH 111 MM, 
MM, M El, M, SI, [H, M, 13 have been investigated 
extensively in QPTs and are recognized to be effective 
and powerful in detecting the critical point. The for- 
mer measures quantum correlations between partitions, 
while the latter measures the distance in quantum state 
space. Therefore, the success of them in characterizing 
QPTs is understood by regarding the universality of the 
critical behaviors itself, that is, the divergent of the cor- 
relation and the dramatic change of the ground state 
structure. Furthermore, as the fidelity depends compu- 
tationally on an arbitrarily small change of the driving 
parameter, Zarnardi et al. suggested the Riemannian 



metric tensor 18 1 , while You et al. suggested the fidelity 
susceptibility [19( , both focus on the leading term of the 
fidelity. In the following, we mainly consider the fidelity 
susceptibility (FS). 

Until now, most efforts have been devoted to the study 
of the global ground state fidelity susceptibility (GFS), 
denoted by \g, which reflects the susceptibility of the 
system in response to the change of certain driving pa- 
rameter. In this work, we study the responses of a sub- 
system, for which we study its FS, the so-called reduced 
fidelity susceptibility (RFS), denoted by Xr- Some spe- 
cial cases have been studied in Refs. [2(| HH, [2?], [28|, 
where the subsystems are only one-body or two-body, 
while in this paper we will study an arbitrary M-body 
subsystem. The motivation for the investigation of RFS 
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is clear in physics. Firstly, it reveals information about 
the change of the inner structure for a system that un- 
dergoes QPT. Secondly, as the existence of interactions 
and correlations, a general quantum system is not the 
simple addition of its different parts, especially in the 
critical region, where the entanglement entropy is diver- 
gent [H, UE3|. Therefore it is significant to investigate 
the behavior of the RFS, as well as the effects of entangle- 
ment on it, in both critical and noncritical regions. And 
our study can be viewed as a connection between the FS 
and the entanglement entropy. 



To study this question, we consider an iV-body Lipkin- 
Meshkov-Glick model (LMG) [II] model, and study the 
RFS for its M-body subsystem. As < \r < Xg 51, 
we consider a more useful quantity, rj = Xr/Xg, an d thus 
7/ e [0,1]. We find that, the behaviors of the RFS, as 
well as 77, are quite different in noncritical and critical 
regions. In noncritical region, the entanglement entropy 
is saturated by a finite upper bound, and the inner cor- 
relation is small, thus the RFS behaves similar with the 
GFS, and the ratio 77 depends on r = M/N but not N . 
However, at the critical point, the entanglement entropy 
tends to be divergent with the increasing of system size, 
and the inner correlations are very strong. Then we find 
the RFS approaches GFS with the increasing of N, and 
j] = 1 in the thermodynamic limit for r 7^ 0. These 
can be understood by considering the divergent of cor- 
relation in second-order QPTs, which is reflected by the 
entanglement entropy. 



This paper is organized as follows. In Sec. II, we intro- 
duce the LMG model and give a brief review of the GFS 
studied in [27]]. Then in Sec. Ill, we derive the RFS in 
the thermodynamic limit and obtain its divergent form in 
the vicinity of the critical point. Then we perform some 
numerical computations, and the results are in perfect 
agreement with our analytical prediction. 
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II. LMG MODEL AND GLOBAL FIDELITY 
SUSCEPTIBILITY 

The LMG model, originally introduced in nuclear 
physics and has found applications in a broad range of 
other topics: statistical mechanics of quantum spin sys- 
tem [3(|, Bose- Einstein condensates |3ll ]. or magnetic 
molecules such as Mni 2 acetate [Hj], as well as quantum 
entanglement l&jl. a nd quantum fidelity [13, [28| . It is an 
exactly solvable [34, l3o| many-body interacting quantum 
system as well as one of the simplest to show a quantum 
transition in the regime of strong coupling. The quan- 
tum phase transition of this model can be described by 
the symmetry broken mechanism, the two phases are as- 
sociated with either collective or single-particle behavior. 
The Hamiltonian of the LMG model reads 

H = -±(S 2 x + 7 S*)-hS z , (1) 
where S a — Y^iLi °a/2 ( a — X ,D, z ) are the collective 



spin operators; a % a are the Pauli matrices; N is the total 
spin number; 7 is the anisotropic parameter. A and h 
are the spin-spin interaction strength and the effective 
external field, respectively. Here, we focus on the ferro- 
magnetic case (A > 0), and without loss of generality, we 
set A = 1 and < 7 < 1. As the spectrum is invari- 
ant under the transformation h —h, we only consider 
h > 0. This system undergoes a second-order QPT at 
h = 1, between a symmetric (polarized, h > 1) phase 
and a broken (collective, h < 1) phase, which is well 
described by a mean-field approach [36]. The classical 
state is fully polarized in the field direction ((c*) = l) 
for ft, > 1, and is twofold degenerate with (V*) = h for 
h<l. 

Before deriving the RFS, we give a brief review of the 
GFS of the LMG model that has been studied in Ref. 
[27|], where the authors employed the Holstein-Primakoff 
transformation and derived the GFS for both phases in 
the thermodynamic limit, 



N 



h 2 (h 2 --/)'' 



Xg (h,j) = < 



.> for < h < 1, 

V(l - h 2 ) (1 - 7) 32(l- 7 ) 2 (l-/ l 2) 2 ' 



(1-7)' 



I 32(ft- 7 ) 2 ^-l) 2 ' 



(2) 



for h>l. 



It has been found that, when h < 1, the GFS increases 
with N and can be viewed as an extensive quantity, how- 
ever, when h > 1 the GFS is saturated with an upper 
bound, i.e. it is intensive. 

III. REDUCED FIDELITY SUSCEPTIBILITY 
A. Thermodynamic limit 

Now we give some basic formulas for fidelity and its 
susceptibility. As the subsystem is represented by a 
mixed state, we introduce the Uhlmann fidelity [421 ] . 

F(p ) p)=tvVp 1 / 2 PP 1 / 2 , (3) 

where p = p (h) and p = p (h + dh) with a certain pa- 
rameter h. If dh tends to zero, the two states are close 
in parameter space, and their Bures distance (4lj is, 

dsl=2[l-F(p,p)}. (4) 

In the basis of p, denoted by {IV^Ij the Bures distance 
can be written as [43| 

1 N J 2 T N / \2 

^ 1 Pn ^ , Pn 1 Pm 

n=l n^tm 



where pi are the eigenvalues of p, N is the dimension of p. 
As FS is the leading term of fidelity, i.e., F = 1 — x^ 2 /2, 
we can get FS for h immediately, 

n—1 ± n.^ni 

(6) 

where dh '■= d/dh. In our study, p and p are just the 
reduced density matrices for ground states. 

In the follows, the A-body LMG is divided into two 
parts, A and B with size M and N — M, respectively. 
We will study the RFS for subsystem A, the reduced 
density matrix is pa- This study would give a connec- 
tion between the RFS and the entanglement entropy [l(| • 
As we know that, the entanglement reflects the correla- 
tion among inner partitions, and our study will reveal 
the effects of these correlations on RFS, especially at the 
critical point. 

Now we introduce the total spin operators for the two 
subsystems, S^' B — J2ieA b a a/^- To describe quan- 
tum fluctuations, it is convenient to use the Holstein- 
Primakoff representation of the spin operators [37j , and 
the first step is to rotate the z axis along the semiclassical 
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magnetization 



Sy I = 

S z 



cos 9o sin 8q 

1 
- sin 6>o cos 9 




(7) 



As presented in [36], 9q = for h > 1 so that S = S, and 
#o = arccos/i for /i < 1. The Holstein-Primakoff repre- 
sentation is then applied to the rotated spin operators 

= M/2 - a^a, 

S A = y/Ma,1yJl-rta/M= (§fj , 

S B = (N — M) /2 - b\ 

S B = VN- Mb 1 yj 1 - fttfe/ (iV - M) = [s B ) f , (8) 

where a (aM and & (6^) are bosonic creation and anni- 
hilation operators for subsystem A and B, respectively, 
and S± B — S^' B ± iSy' B . After this transformation, 
the LMG Hamiltonian is mapped onto a system of two 
interacting bosonic modes a and b. For fixed r = M/N, 
the Hamiltonian can be expanded in 1/N. Up to the or- 
der (1/A)°, one gets H = NH^ +H®> +0 (1/N) with 



1 — 2h)/4, where to = cos6*o, and 



i/ (o ) = _l+7 +A t VA T^ 1 - 

4 

where A = (a, 6) , and 

v _2/im + 2 - 3m 2 - 7 



AtWfAtl 



h.c. 



(9) 



W 



7 — to 



T V 7 " (1 - T) 

\/t(\-t) l-T 



(10) 



where I is a 2 x 2 identity matrix; to = h in broken phase 
and to = 1 in symmetric phase. The bosonic Hamiltonian 
can be diagonalized by Bogoliubov transformation and is 
useful in deriving the reduced density matrix. As shown 
in [H, H^, H(| , the reduced density matrix for eigenstates 
of a quadratic form can always be written as pa = e 
with 



H 



U = k + K\a) a + k 2 (a t2 + a 2 ) . (11) 

/tj (i = 0, 1, 2) can be determined by using (loj 

trp^ = 1, tr (p^a^a) = (a^a) and tr (p^a* 2 ) = (a^ 2 ) ■ 

(12) 

where (CI) = (ip g \fl\i(jg), \ip g ) is the ground state, Then we 
can diagonalize pa by Bogoliubov transformation. How- 
ever, in this paper we will adopt another method to di- 
agonalize pa, as shown in Ref. [Ill ], pa is written in the 
bosonic coherent state representation 



\pa\$) = if exp 



G 4 



1 



G~> 



x exp 



(^* + 00 
•-00 



where a\<f>) = 4>\<j>)\ K = y/(l + G++) (1 - G ) is de- 
termined by the normalization of pa', G ++ and G are 
Green's functions defined as 



C++ = ((at 
G"- = ((at 



(13) 



Then p" 4 can be diagonalized by the following Bogoliubov 
transformation, 



g = cosh tpa + sinh <pa^ 

P + Q , P-Q t 

= a H a 1 

2 2 



(14) 



with PQ = 1, PG++ = pQ, and QG = -pP. The 
Green's functions can be obtained by diagonalizing the 
bosonic represented Hamiltonian Q, 



G++ = l + (l/a-l)r, 



(1 



where 




a)T - 



for h>l. 



(15) 



(16) 



for < h< 1. 



The diagonalized p A reads 



1 



(17) 



where the pseudoenergy e = In [(/1 + 1) / (jj, — 1)] with 
H = a" 1 / 2 ^™ + (1 - t)] [r + a (1 - t)]. 

Now we can derive the RFS, of which the first term 
involves only the eigenvalues of pa, and the second term 
involves both the eigenvalues and the eigenvectors. The 
eigenvectors of pa is the number state \n): g^g\n) — n\n), 



and the term \(tp n \dhipn 
lated by using 



\(n\dhm)\ can be calcu- 
(18) 



1/ iq \|2 \{ n \dh.g^ g\m)\ 

\(n\d h m)\ =^— -3^-. 

(m — nj 

Then we write the RFS explicitly, 

^ ( ^ 7 ' r) = 4T^-T) + ^Tr + ^ ( ^ ' 

(19) 

where ip = arctanh [(p - G++) / (p + G++)], O = 
arccosft for ft < 1 and $o = for h > 1. Thus the 
last term of the above expression only takes effect in the 
broken phase. We emphasize that, in the broken phase 
h < 1, we should perform a rotation ([7]) at first. 



4 



8000 
7000 
6000 


N=512 o 
- V=1024 • t, 
\ 2048 o 


_ 5000 
eg 




7 4000 


f i 


N 3000 


♦ <> 

&\ : 


2000 


J5 * : 


1000 






0.8 0.85 0.9 0.95 1 1.05 1.1 



FIG. 1: RFS as a function of /i at 7 = 1/2 and r = 1/2. 
The peaks approach the critical point and become sharper 
and sharper with the increasing of N. 



We can express it farther as 

Nt 



Xr (h,J,r) = 




4G++ (1 - h 2 ) 



for < ft< 1, 
for ft > 1, 



where 



d h In - 



G++ 



(20) 
(21) 



4(m 2 -1) 4(/i 2 + l) 
In the vicinity of the critical point, the RFS diverges as 

Xr/Noz (1 - /J) -172 , for < h< 1, (22) 
Xr oc (1 - ft)~ 2 , for ft > 1, (23) 

and this is the same with \g- Additionally, we show the 
ent ang lement entropy £ = — tr(plnp) that was derived 
in dflHII, 



£ = 



iln^- 



1 



iln^- 



1 



;ln2. 



(24) 



where x — 1 when h < 1 and x = when h > 1, 
the In 2 term comes from the two-fold degeneracy of the 
ground state in the broken phase, and this degeneracy 
is lifted for finite N. The entanglement entropy diverges 
as (1/4) In \h — 1| around the critical point, and is nearly 
independent with TV in noncritical region. 



B. Finite size cases 
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FIG. 2: A comparison between r\ (a) and £ (b) as a function 
of ft at 7 = 1/2, t = 1/2 for various system sizes. At the 
critical point, r\ tends to 1 while £ is divergent. 



\*9 



N 

^2 C m \J, -J + to), 



(25) 



where C m is the coefficient to be determined numer- 
ically. We hope to write \J,—J + to) in the form 
of I Ja, TOa)| Jb) where Ja = M/2 and Jb = 

(N — M) /2 correspond to the two local systems. Since 
\J,-J + m) = ^{2J-m)\/ (2J)!m! (S+) m \J, - J), and 
the ladder operator S+ = S+ + S 13 . Then the ground 
state is 



\*l>9 



where 



N 2J A 

= ^2J^C m y/H(p; 2J,2J A ,m)\JA,-JA+p) 

m— p—0 

®\Jb,-Jb +m-p) (26) 



To perform numerical computations, we should derive 
the reduced density matrix for pa in finite size case. The 
LMG model is of high symmetry in interaction, and the 
ground state which is the superposition of the Dick states 
lies in the J — N/2 section 



H(p;2j,2j 1 ,m)= K P ^ (27) 

is the so called Hypergeometric distribution function. 
And the matrix element of pA is 
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N 

(pA) p<q = J! CmC* +m _ pV /H (P: 2 J: 2 ^A, TTl) 

x v/H (g;2J,2J A ,g + m-p). (28) 

By using the exact diagonalization method, the RFS 
as a function of h for fixed r is computed and shown in 
Fig. fTJ). As one can see that, the peaks of the RFS ap- 
proach the critical point and become sharper and sharper 
with the increasing of N. The RFS in the symmetric 
phase (h > 1) has an upper bound, however, in the bro- 
ken phase (h < 1) the RFS increases with the total spin 
number iV. Thus we address that, the RFS is extensive 
in the broken phase, in which the LMG model is of collec- 
tive behavior, while is intensive in the symmetric phase, 
in which the LMG model behaves like a single particle. 
This is similar with the GFS [27j]. 

As < Xr < Xgi we "will focus on a more useful quan- 
tity r)(r, h) ee Xr (h, 7, r) /xg (h, 7) and study its prop- 
erties in critical and noncritical regions. With Eqs. 
(|19p . we find that in the thermodynamic limit 

lim r] (t, h) = 1, (29) 
h— *1 

for any non-vanishing r. To verify our prediction, we 
show the analytical and numerical results in Fig. 
As one can see that, at the critical point, the RFS ap- 
proaches the global one, i.e. r] tends to 1, and at the 
same time, the entanglement entropy, i.e. the inner cor- 
relation between subsystems A and J3, is divergent with 
the increasing of N. When h is away from the critical 
region, the inner correlation decreases dramatically, and 
then rj depends on r but not the total system size N as 
shown in Fig. ©. 

As demonstrated in Ref. [28j, when there are no cor- 
relations between partitions of a system, for example an 
TV-body system represented by a product state that reads 

JV 

IV>W) = (g)l&W>, (30) 

i=l 

if we denote a one-body reduced fidelity as F r , the rela- 
tion between the global and the reduced fidelities is 

n 

F g (h,S) = l[F;(h,S). (31) 

1=1 

and thus we have Xg — X)-j=i X r i moreover, if the system 
is of translation symmetry, we have Xg = ^Xr- If there 
is entanglement between partitions, we have no such re- 
sults, especially in the critical point, the entanglement 
is divergent, and then Xg/Xr = 1 in the thermodynamic 
limit. This is some kind of effect of the inner correlations 



on the susceptibility of the system states. However, we 
address that our results are based on a high-dimension 
model, actually there are interactions between any two 
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FIG. 3: rj as a function of r with 7 = 1/2, at h — 0.6 (a), 0.9 
(b), 1.0 (c) and 1.1 (d). We see that rj is nearly independent 
of N when h is away form the critical region. 



particles in the LMG model. We think it is deserved to 
study the RFS for a contiguous block in a low-dimension 
model, for example, the XY model in which the interac- 
tion is just between neighboring sites. Thus the correla- 
tion between a block and its complementary part takes 
effect only on the boundary, and the results for rj maybe 
different. 



IV. CONCLUSION 

In conclusion, we derive the RFS analytically in the 
thermodynamic limit for a fixed r. To analyze the effects 
of the inner correlations on the RFS, we study the ratio 
rj = Xr/Xg combined with the entanglement entropy in 
both critical and noncritical regions. Our results give a 
clear picture for understanding the effects of correlations 
on the response. In the critical region, with the increasing 
of N, the entanglement entropy tends to be divergent 
and 77 approaches 1, while in the thermodynamic limit, 
77 ee 1 for t ^ 0. This indicates that, the sensitivity of 
the subsystem is equal to the global one. In noncritical 
region, the RFS behaves similarly with the GFS, and ij 
depends on r but not N. 
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